Spatial Coherence of a Polariton Condensate 
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We perform a Young's double-slit experiment to study the spatial coherence properties of a two- 
dimensional dynamic condensate of semiconductor microcavity polaritons. The coherence length 
of the system is measured as a function of the pump rate, which confirms a spontaneous build-up 
of macroscopic coherence in the condensed phase. An independent measurement reveals that the 
position and momentum uncertainty product of the condensate is close to the Heisenberg limit. An 
experimental realization of such a minimum uncertainty wavepacket of the polariton condensate 
opens a door to coherent matter-wave phenomena such as Josephson oscillation, superfluidity, and 
solitons in solid state condensate systems. 

PACS numbers: 



Simple, yet profoundly connected to the foundation of 
quantum physics, the Young's double-slit experiment has 
been a benchmark demonstration of macroscopic spatial 
coherence - off-diagonal long range order (ODLRO) of 
a macroscopic number of particles [l| - in Bose-Einstein 
Condensation (BEC) of cold atoms [2, & H . Recently, a 
similar phase transition has been reported for the lower 
branch of exciton-polaritons (LPs) in planar semicon- 
ductor microcavities 0, H, 0, S, H \m, EL El, and 
supporting; theoretical frameworks have been developed 
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Interestingly, LPs are 
free particles in a two dimensional (2D) system where 
genuine BEC exists only at zero temperature in the ther- 
modynamic limit 2l|, |22j . A quasi-BEC can be defined 




for a 2D system of a finite size if a macroscopic number of 
particles occupy a single ground state and if an ODLRO 
is established throughout the system 23j, [24J . Yet in the 
LP experiments to date, the system size is ambiguously 
defined by the spot size of the pump laser, and there 
is no quantitative study of the relation between the size 
and the coherence length of a condensate [Hi 11 1. In 



this work, we perform a Young's double slit experiment 
on a LP gas to measure its spatial coherence properties 
across the phase transition, and compare the measured 
coherence length with the condensate size. We also mea- 
sure the position-momentum uncertainty product of the 
condensate and compare it to the Heisenberg limit. 

A sketch of the setup is shown in Fig. 1. The mi- 
crocavity sample is first magnified by a factor of 37.5 
and imaged to a plane A, which is in turn imaged by a 
lens II to a charge-coupled-device (CCD) at plane C for 
measurement of spatial distribution. For the double-slit 
experiment, we insert a pair of rectangular slits at plane 
A, and move the lens II such that the image of plane 
A (denoted by plane B) is a distance D behind plane 
C. Effectively, we observe on the CCD the interference 
pattern of the LP emission passing through the double- 



FIG. 1: A sketch of the double-slit experiment setup. Upper: 
the LP spatial distribution is imaged to a CCD at plane C. 
Lower: a double-slit is inserted at plane A and imaged by lens 
II to a virtual plane B. The CCD at plane C, a distance D 
from plane B, captures the double-slit interference pattern. 



slit. In our experiment, D — 6.7 cm, the width of the 
slit image at plane B is 8 = 53 /im, and the average 
wavelength of the LP emission is A ~ 778.5 nm. Corre- 
spondingly, the Fresnel number = 0.05 <C 1, thus the 
far-field condition is satisfied at plane C. When mapped 
onto the sample surface, the slit width seen by the LPs is 
Ar w 0.5 ^m, which is less than the intrinsic coherence 
length £o ~ 1 ^ni of a single LP [25| , and much less than 
the LP system size of 5 — 10/im. Hence neglecting the 
variation in LP distribution within each slit, we obtain 
the intensity distribution on the CCD camera [261 ] : 



I(x) = h(x) + I 2 (x) + 
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The subscript i = 1, 2 denotes the slit number, r, are the 
x-coordinates of the slits on plane B (Fig. [1]), a; is the 
x-coordinate on plane C, Xo is the center of the double- 
slit on plane C, Efc) is the LP field amplitude at slit 
i, ktot is the free space average wavenumber of the LP 
emission, and d is the separation between the images of 
the two slits at plane B. Ii is the intensity distribution 
if only slit i is open. <^i2 is a fixed phase difference of 
LPs between the two slits. <f>{x) is a varying phase close 
to the path length difference from the two slits, giving 
rise to a cosine modulation on the far-field intensity dis- 
tribution. After proper normalization, the amplitude of 
the cosine modulation equals the first order coherence 
function g^'(r = |ri — r%\). Using six sets of double-slit 
with varying slit separations r, we measured (r) from 
1.3 /xm, close to the intrinsic coherence length of a single 
LP, up to 8 /im, close to the LP system size. By varying 
the pumping intensity, we studied the characteristics of 
g^\r) across the phase transition. 

The sample we investigated has a A/2 GaAs cavity 
sandwiched between Gao.865Alo.i3sAs/AlAs distributed 
Bragg reflectors. Three stacks of quantum wells (QWs) 
are placed at the central three antinodes of the micro- 
cavity, each stack consisting of four 6.8 nm-wide GaAs 
QWs separated by 2.7 nm-wide AlAs barriers. We pump 
the sample with linearly polarized pico-second mode- 
locked Ti-Sapphire laser. At an incidence angle of 50° 
from the sample growth direction, the laser is resonant 
with the exciton-like LP modes. The sample is kept at 
Tiattice ~ 4 K. The cavity-photon energy at zero in-plane 
wavenumber is ~ 7 meV above the bare QW exciton reso- 
nance, corresponding to an optimal detuning for thermal- 
equilibrium condensation of the LPs 1121. The threshold 
pumping density is Pth ~ 300 W/cm 2 [271 ]. 

A typical interference pattern observed at a pump rate 
above the condensation threshold is shown in Fig. [2j\. 
Distinct interference fringes are readily observed imposed 
on a sine function distribution. To obtain g^{r), we in- 
tegrate over a narrow strip along the y-axis and fit it with 
Eq.[TJ <?W, 0i2, and I-E2I are free fitting parameters, 
while xq, d, X and X c are estimated from experimental 
parameters with a 10% allowed variation. As shown in 
Fig. [2B, Eq. ED fits the data very well for P > P th . At 
P < P t h, the interference patterns are barely observable 
or non-existing; one example is given in Fig. 

In Fig.[3]A_ we show the increase of g^'(r) with P/Pth 
for a few slit separations r. A salient feature in Fig. 
is that, there is a jump in g^'{r) when the pump rate 
is increased above a condensation threshold, even at r 
up to the laser pump spot size of ~ 8 /im. When the 
pump spot size is increased to ~ 20 /im, g^ x \r) > 0.3 
was also observed up to r = 20 /im, limited again by 
the pump spot size [28(. This demonstrates the sudden 
appearance of macroscopic coherence above the conden- 
sation threshold. Another feature is that, the increase of 
g^'(r) with P/Pth is slower at larger r. This shows that 
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FIG. 2: (A) Raw image of the interference pattern, slit sep- 
aration r = 2.7 /im, P/Pth = 7. (B) Measured (symbols) 
intensity I(x) for r = 2.7 /im, Pj Pth=6-7 and fitting by 
Eq. Q] Fitted g {1) (r) = 0.560 ± 0.006. (C) Same as (B), 
for r = 2.7 /im, P/Pth = 0.5 and g w (r) = 0.09 ± 0.02. 



the macroscopic coherence is built up gradually through- 
out the pump beam spot size when the phase space den- 
sity of LPs is increased. 

To study the spatial coherence properties quantita- 
tively, we plot in Fig. [3)3 how g^(r) decays with r, 
and define a coherence length r c as g^(r c ) = 1/e. 
The normalized pump power dependence of r c is shown 
in Fig. [3p. As a reference, for a classical Maxwell- 
Boltzmann (MB) gas in thermal equilibrium: 



g (1) (r) 
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Here Ap is the thermal de Broglie wavelength, m is the 
mass of the particles and T is temperature. For the cur- 
rent system at T — 4 K, A t /^/tt w 1.9 /im. 

Below condensation threshold density, we measured 
a r c < At / \/tt, since the system is far from ther- 
mal equilibrium. In fact, the finite g {1) {r) - 0.15 at 
r = 1.3 /im is consistent with the intrinsic coherence 
length £0 ~ 1 M m °f a single LP due to its finite lifetime 
[25j ]: <7o ( r ) = ex P(~p') = 0.19. Above condensation 
threshold, we found r c 3> At/v 7 ^"- Fitting of the data 
with Eq. [2] with temperature T as a free parameter also 
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FIG. 3: (A) g^'(r) vs. P/Pth for different slit separations r 
as labeled in the figure. (B) g^'(r) vs. r for different pump 
rates P/Pth as given in the legend. The symbols are measured 
g^'(r). The solid lines are fittings by Eq. [3j The dashed line 
at P/Pth — 7 is a fitting by Eq.[2] The dash-dotted line marks 
where g (r) = 1/e. (C) r c , the 1/e decay length of g (r), 
vs. P/Pth (circles) and Afc vs. P/Pth (triangles). 



fails to describe the data (dashed lines in Fig. 03). 

Since g^'(r) is the Fourier transform of the momen- 
tum distribution /(fc), we resort to the actual momen- 
tum distribution of the system. It was found that above 
the condensation threshold, the LPs become highly de- 
generate in the states with the lowest kinetic energies 
(e.g., Fig. 2 in Ref. [3). Their momentum distribution 
deviates from the MB distribution, but follows well the 
Bose-Einstein (BE) distribution with chemical potential 
I A* I *C k B T. In this quantum degenerate limit, f(E ~ 
k B T) = [exp( j^) - l]" 1 « (e - l)" 1 « /(0) « k B T / ^ 
most of the emission comes from LPs with E ~ 0. Hence 
we can obtain the following approximate form of g^(r): 

gW(r) oc ^ 2D \f(k)g{k))^H{f{k)) 
k B T 



E{k) - f i' 



xKo{r VEfEEr y (3 ) 

At 



Here k is the LP's in-plane wavenumber, g(k) is the 
constant momentum density of state, T^ 2D \f{k)) and 
H(f(k)) denote the 2D Fourier and Hankel transform 
of f(k), respectively. As shown in Fig. [3)3 (solid lines), 



Kq (x) , the modified Bessel function of the first kind, fits 
very well the measured g^'ir) for P/Pth > 1- 

If we extrapolate the e _Tr decay length of (r) from 
the fitting and compare it to At at 4 K, it is 2.5At at 
P = l.2P th , and 7.9A T at P = 6.7P th . 

In Fig. we compare the pump rate dependence of 
r c and Afc, the 1/e width of the measured first-order 
coherence function g^\r) and the measured momentum 
distribution function /(fc), respectively. At pump rates 
lower than the condensation threshold, r c is ~ 1 /an, 
limited by the intrinsic coherence length £o of a single LP, 
while Afc is ~ 2 /im -1 due to the slow energy relaxation 
dynamics of the LPs. The product r c ■ Afc is close to 2, 
the value expected for a thermal MB distribution. When 
P increases toward Pth, more injected LPs relax to the 
lower energy states and Afc gradually narrows, but r c 
is still limited by £o> hence a decrease in r c ■ Afc. Once 
above the threshold, there is a sudden increase of r c by 
more than five fold up to ~ 6 /im, which manifests the 
spontaneous build-up of a global phase among the LPs 
due to the phase coherent stimulated scattering of LPs 
into the ground state. Correspondingly, Afc is reduced by 
about four fold since the LPs form a quantum degenerate 
Bose gas. Further increasing the pump rates, r c decrease 
slightly while the momentum distribution is broadened, 
potentially because stronger LP-LP scattering at high 
densities introduces condensate dephasing [19| , 

Finally, it is instructive to compare Afc with the mea- 
sured condensate size 0. We have consistently observed 
an abnormally slow increase of the condensate size in 
comparison to the spot size of a photon laser based on 
electron- hole pairs (Fig. 6C in Ref. 0). 

Due to the discrete jump in quantum efficiency at the 
condensation threshold (Fig. 1 in Ref. Q), the emission 
in a condensate region is much brighter. Since the pump 
beam has a Gaussian spatial profile, the center of the spot 
reaches a threshold first, leading to a sharp decrease of 
the emission spot size at Pth, for both a LP condensa- 
tion and a photon laser. At P > P t h, the emission spot 
size measures the area which reaches the threshold. In a 
photon laser, threshold density is determined by the local 
density of electron-hole pairs independent of the system 
size. Hence the 1/e spot size u) can be estimated as: 



w(P/Ptfc)=WpWl-l0g a (l + -^) 



Pth 
p 



(4) 



where lo p is the pump spot size. Eq. [5] describes very well 
the photon laser data (Fig. 6C in Ref. [9(), but fails to ex- 
plain the data of a LP condensate (Ref. [9] and Fig.fiJ'V). 
Here we propose that for a LP condensate, uj (P/Pth) re- 
flects the size of the condensate in which a condensation 
threshold is satisfied. Then the critical LP density n c (u>) 
increases with the system size u>, and the pumping rate 
P/ Pth needs to be modified as „ ("(p/p fh )) ; where u> c 
is the condensate size at P = P t h- As a simplified model, 
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FIG. 4: (A) System size 2to vs. pump rates P/Pth- Symbols 
are the data. The solid line is a fitting by Eq. [5j with uj p = 
13.3 ± 0.1 fim, uj c = 6.1 ± 0.1 /im. The dotted line is a 
fitting by Eq. U for comparison, with uj p — 11.9 ± 0.2 /xm. 
(B) The position and momentum uncertainty product Ok ■ 
o r vs. P/Pth- The dash-dotted line indicates the minimum 
uncertainly of at ■ ay = 1/2. 



consider tt^ for a 2D boson gas confined in a finite size 



L = 2uj 



w{P/P t 



n c {uj) 



\n{j^-), then we obtain: 



th 
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P th ln(2^/A T ) , 
P ln(2w c /A T )' 
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Eq.[5]fits the data very well (solid line in Fig.[UA), with uj p 
and lo c as fitting parameters. This suggests that 2tu mea- 
sures the size of the coherent condensate above threshold. 

With the distribution functions in both spatial and mo- 
mentum domains, we can evaluate how well the system 
can be described by a single-particle wavefunction. The 
standard deviation and a T are calculated from the mo- 
mentum and spatial distribution data f{k) and f(r) at 
P > Pth, respectively; their product is compared to the 
Heisenberg minimum uncertainty limit in Fig. The 
sharp decrease of • oy at P ~ P t h indicates that a large 
number of the LPs in the system condense into a sin- 
gle quantum state. Deviation from the Heisenberg limit 
shows that there are some thermal LPs coexisting with 
the coherent condensate. The slight increase of Cfc • ay at 
P/Pth > 1 may be caused by condensate depletion due 



3 SI 



to LP-LP interactions at high densities 

In conclusion, we studied the spatial coherence of a 
microcavity polariton condensate. A Young's double-slit 
setup is implemented to measure the first order coher- 
ence function g^(r) of the LPs. The system acquires 
macroscopic coherence above a condensation threshold, 
manifested as a sudden jump of g^'(r). The observed 



g( x \r) vs. r is well described by the Fourier transform of 
a degenerate Bose-Einstein distribution in the momen- 
tum space. The coherent condensate expands from the 
central region of the pump spot to the full pump spot size, 
and a slow growth of the condensate size is well under- 
stood by a simple model of quasi-BEC with a finite size. 
We also confirm the position and momentum uncertainty 
product cTfcoy of the LPs decreases toward the Heisenberg 
uncertainty limit above a condensation threshold. 
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